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ABSTRACT 



In order to provide a theoretical tool well suited for use in 
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characterizing the stability margins (e. g. , gain and phase margins) of 


U 9 N multiloop feedback systems, multiloop input-output stability results 
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generalizing the circle stability criterion are considered. Generalized 


N 

O 


K z ca 

O M QS 

MU o 

H OS 

■< •< U 

N K — .U 

M t Vi 

“4 x: u 
<« H o 
as M o 

W 1-1 f-i 


a: a 
u ^ 

H 
a. VI 

Sc.- 

mi C 

M a. 


w 

e o *4 


... a] 

I" 

M 

... OB 

U 

u 

rv u 

<N 

•“ cs. 

I M 
« U 

u 

I u 
< s 

tr, H 


X* 

0) 

•3 M 

x: e 
V 

f'i a. 
oj s- 
« N. 
rc m 
«— o 

vj 

M O 
W5 X 
M 

-J O. 


z a. z rst 
— o <« <-> 


conic seccors with "centers" and "radii" determined by linear 

dynamical operators are employed to enable an engineer to specify the 

stability margins which he desires as a frequency-dependent convex set of 

modeling errors- - including nonlinearities, gain variations and phase 

variations- - which the system must oe able to tolerate in each feedback loop 

without instability. The resulting stability criterion gives sufficient condit • 

ions for closed- loop stability in the presence of such frequency- dependent 

modeling errort, even when the modeling errors occur simultaneously in 

all loops; so, for example, stability is assured as loop gains and phases 

vary throughout a "set of non- zero measure" whose boundaries are 

frequency-dependent. The stability conditions yield an easily interpreted 

scalar measure of the amount by which a multiloop system exceeds, or 
tails snort of, its stability margin specifications. 
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I. INTRODUCTION 


A key step in the synthesis of robustly stable feedback systems is 
the characterization of a set of feedback laws that are stabilizing for every 
element of the set of possible plant dypjimics. This type of information is 
precisely what'is provided for single -loop feedback systems by such 
input-output stability criteria as the Nyquist, Popov, and circle theorems. 
Indeed, the practical merit of classical feedback design procedures in- 
volving Nyquist loci, Bode plots, and Nichols charts is in a large measure 
directly attributable to the fact that these design procedures provide the 
designer with sasily f.nterpretabie characterizations of such sets of ro- 
bustly stable feedback laws. For single-loop feedback systems, these 
stability theorems enable engineers to meaningfully characterize the 
tolerable amount of gain and phase variation in the loop at each frequency, 
and even the tolerable amount of unmodeled nonlinearity. These tolerances 
of modeling error are in broad terms what we call stability margin , 
classical gain and phase margin being two familiar measures of stability 
margin for single feedback loops. Although multiloop generalizations of 
the Nyquist stability criterion have been developed (e. g, , [ IJ -[3]), it 
has been difficult to meaningfully relate tho conditions of these multiloop 
criteria to tolerance of open-loop modeling error except in special cases 
such as diagonally dominant systems, normal systems, and systems in 
which feedback loop gains vary only over certain "sets of zero measure" 

[ 4] , The results of the present paper are intended to address the need 
for an improved method for characte rizing the stability margins of multi - 
loop feedback systems. 

In broad and imprecise terms what seems to be necessary to meaning 
fully characterize multiioop stability margins is a stability criterion that 
guarantees stability for every multiloop feedback operator within a given 
"frequency-dependent ball" in an appropriate space of input -output relations, 


this ball being centered at the system's nominal "open -loop gain" operator. 

It is important that the size of this ball be permitted to be frequency- 
dependent so that one can account for frequency -dependent variations 
in the precision of mathematical models such as result from such ubiquitous 
effects as singular perturbations* hysteresis, imprecisely known time- 
delays, or any sort of unmodeled dynamics. Also, since in general 
one may expect modeling imprecision in certain feedback loops to be large 
relative to other loops, it should be possible to specify that this ball be 
somewhat egg-shaped, having different diameters in the various "directions" 
corresponding to the "gains" of individual feedback loops. So, perhaps 
the necessary ball of stable multiloop feedback operators could be better 
described as a "frequency-dependent egg. " All of this is of course too 
vague and imprecise to be of immediate use«-what is needed is a stability 
criterion dealing with a precise mathematical description of this frequency- 
dependent ball (or, "egg") and of the space of operators in which is is 
embedded, 

Stability results in this general spirit are provided by the Zames- 
Sandberg input-output stability theory [ 5] -[ 6] , Sandberg's frequency- 
domain stalility criterion[ 6 ] for systems with multiple nonlinearities can 
be interpreted as guaranteeing stability for a collection of nonlinear feed- 
back operators inside a spherical (i. e. , not "egg-shaped") non-frequency- 
dependent ball centered at the identity operator times a scalar; the now 
well-known circle stability criterion emerged in [7] as a special case of 
this result. Zames' conic sector stability theorem [ 5, Theorem Z] is an 
abstract generalisation of Sandberg' s criterion that makes the connection 
with balls of stable muitiloop feedback operators even more transparent; 
the conditions of Zames' theorem involve conic sectors which it happens 
are simply spherical balls, centered at the identity operator times a 
scalar in an extended normed space of input-output relations. 
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In tho past fifteen years frequency domain stability criteria based 
on the Zames -Sandberg theory have been improved and refined in many 
significant way?. Reference [8] provides a good overview of much of 
this work, Reference [9] develops similar results in a Lyapunov setting, 
References [lO] - [12] and the additional references cited therein describe 
many stability results developed specifically for interconnected (i, e, , 
multiloop) nonlinear .‘lystems. However, the previous literature in this 
area has focused primarily on nonlinear stability; though Zames [5] makes 
some key suggestive remarks about the broader implications conic sector 
results regarding imprecisely modeled systems. The stability margin 
implications of the results have not bean stressed and no results based 
on the Zames -Sandberg theory have been published which address the 
need for a frequency-dependent characterization of multiloop stability 
margins. 

The main objective of the present paper is to present a multiloop 
input -output stability criterion that is tailored to the task of multiloop 
feedback stability margin analysis. Our main result ( Theorem 1) shows 
that multiloop stability margins --including tolerance of unmodeled non- 
linearity and of dynamical modeling errors of frequency-dependent 
magnitude --can be directly related to open-loop system frequency -respon se 
quantities. The results, expressed in terms of the ’’singular values" of 

* 

certain matrices, are observed to yield an easily interpretable scalar 

measure of a systems "excess stability margin" --i, e. , of the amount by 
which a multiloop feedback design exceeds its stability margin specifications. 
A related result described in [21] makes use of Theorem 1 in. generating 
generalized conic sector bounds for characterizing the sensitivity of 
multiloop systems to large dynamical modeling errors of frequency- 
dependent magnitudes; in effect the result of [21] provides a nonlinear 
multiloop generalization of classical M-circle ideas. 
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The role of singular values in connection with stability was first 
noted by Sandberg [j'l , though Sandberg does not specifically use the term 
singular value. Earlier versions of the results in present paper and thoir 
connections with stability margin analysis were first reported in [lE] -[14]. 
Stressing the use of singular values, Doyle [15] establishes important 
geometric connections between multivariable Nyquist criteria and the 
stability margin results of [ 13], [ 14], and the present paper; additionally, 
Doyle [is] and Stein and Doyle [16] cite a number of illustrative examples 
that present a compelling case for the use of results of this type in the 
analysis of multivariable feedback stability margins. MacFarlane and 
Scott-Jones [25] discuss at length the relationships between the eigenvalues 
and the "principal gains" (i. e. , singular values) cf a multiloop-system's 
transfer matrix, Nuzman and Sandell [23] establish some inequalities 
relating these results to the singular values of the return difference matrix 
for multiioop systems and discuss the connection with the guaranteed stability 
margins of full-state feedback linear optimal regulators (see [24] ). The 
paper by Sandell [l7] discusses in broad and simple terms the role of 
singular values in coping with modeling imprecision in a wide range of 
engineering ?.nd numerical problems. Numerical aspects of singular value 
computation axe surveyed in [ IS], wherein sophisticated and widely avail- 
able computer routines for singular value computation are aleo .referenced. 

n. NOTATION 

T T 

The following notation is used; A and x denote respectively the 

transpose of the matrix A and the vector x; A and x"* denote the complex 

T T 

conjugate of the matrix A and the ve.ctor x respectively; the determinant 
of a matrix A is denoted det(A); the Euclidian norm, of a vector x is 
!]x|l jg = '/ X X ; R^denotes nonnegative real numbers; the functional 
norm jixjjT- and inner product < x^, -’^2 defined for functions x: 
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U) 


R . - R*' as 

T 


lu II T = ^ 


where for any Xj^ and X2 


< Xj^, ^2 " I • 


(2) 


The space L2^(R^, R^^) is defined as 


= {*= II ^ IIt ^" ''’■''^^.1 


( 3 ) 


where the symbol V means "for all, " We define I-<2 q = U 

n 

Laplace transforms are denoted by capital letters, e. g, , X(s) denotes the 
Laplace transform of x(t). 

Given any matrix A, the square-roots of the eigenvalues of A'A are 
called the singular values of A, For any matrix A, we use the notation 
a (A) to denote the largest singular value of A and a . (A) to denote the 
smallest singular value of A. Singular values are always nonnneg'ative real 
numbers since A A is always positive semidefinite, 

A functional relation is a mapping of functions into sets of functions; 
for example, a dynamical system mapping inputs in L2^ into outputs 
in L_ defines a relatjon (e, g. , [ 5 ]). An operator is a special type of 
relation which maps each input function into exactly one output function, 
i, e, , into a set with exactly one element. All functional relations consi- 
dered in this paper are implicitly assumed to be mappings of L2^ into lj2g* 
A relation H is said to be nonanticipative if for all t^, the output (Hx)(tQ) 
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does not depend on x(t) for t > tQ. We say that a relation H is Lq^- stable ^ 

if there exists a constant k < “ such that for all x c Li^ and all t e R 

2e f 

It K X !1 ^ < k jl X 11^ . (4) 

Generalizing some of the conic sector conditions of Zam«s [5], 
we employ the following definitions which are a special case of the 
generalized sector conditions of [l3], [19], [20]. Given an operator H, if 
there exist operators C, R, and S such that 


l|S(y-Cx)ll^ < llRxIlJ - £(11x11^ + llyll^) (5) 

for all y = Hx, all x, all T, and some e > 0, then we say "H is strictly 
inside the conic sector with center C ind radius (R, S)"; equivalently, 
v/e write 

H strictly inside 1.2g-Cone(C, R, S). ^ (6) 


Given a relation G , if 


ll5(y- > l|R== ll? 


(7) 


^For nonanticipative operators, stability as defined here is equivalent 

to the usual notion of L 2 - stability, e, g. , [8j. 


2 

The "strictly inside 
to be a special case 
in [ 5 ]. 


" conic sector condition of [5] can be demonstrated 
of (6); however, the term ellyjl^ is only implicit 
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for all X = «Gy, all y, and all t, then we say "the inverse relation yf -G 
is outside the conic sector with center C and radius (R, S)"; 
equivalently, we" write 



puts ide 


^2e" S) . 


( 8 ) 


The notation col(x^, . . . denotes the column vector 


r“ X, -1 


col • ■ • » 




(9) 


The relation diag(Hj^, • • > > is defined by 


diag (Hj^, . . » , l^)*c ol(Xj^ , • • • , ~ c ol (H , < i • , . (10) 


3 I 

Following [13], [l9], [20], the notation (-G) is used for the inverse of the 

relation -G ; i. e. , (-G)^ is the relation which maps each y? into 

the set of functions xe such that y = - G x. The inverse relation (-G)^ 

always exists even for operators G for which the inverse operator, 

denoted (-G) ^ , ‘does not exist. 
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HL PROBLEM FORMULATION 


Our results concern the input-output stability of systems consisting 
of a dynamical linear time- invariant (IjTI) interconnection of N imprecisely 
modeled components, including imprecisely model^td LTI components and 
nonlinear time -varying (NTV) components as well as dynamical nonlinear 
components comprised of interconnections of LTI and NTV subcomponents. 
The system equations thus take the following form (see Fig. 1); 

components 


^i 


= H. X 

^ i 1 


(i = 1,...,N) 


( 11 ) 


dynamical LTI interconnection 


X(s) = -G(s) (Y(s) + V(s)) + U(s) 


( 12 ) 


where 


Y(s) = col (Yj(s),...,Yj^(s)) 
X(s) 

U(s) 

Y(S) 


col (Xj(s),...,Xj^(s)) 
col (^^^(s), . t . , Ujg(s}) 
col (Vj(s),...,Vj^(s)) 


(13) 

(14) 

(15) 

( 16 ) 


e(s) = 


Gj^(s) 


Gj^l(s) 


« « • 
• • • 


Gin(s) 




(17) 
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[COMTOI^Ej^S 



INTERCONNECTION 










Fig. 1 The System 
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The endogenous variables yj^{t) and Xj(t) are the system ''outputs" 
and the exogenous variables u^(t) and Vj^(t) are the system "inputs". Each 
of the "compbhents" H. may itself be a multi,- Input-multi-output a (MIMO) 
system In general, though our results arc most easily used and inperpreted 
when the components are slngle-lnput-slngle-output (SISO). 

We assume that for each of the imprecisely modeled components 
H. we have a crude approximate LTI model C . and that LTI operators R . 
and can be found such that the modeling error in each is bounded by 
a generalized conic sector condition such as (6); i. e, , 

H. strictly inside - Cone{C., R., S.) (18) 

n# X wC MP X M X X 

for is 1, , . , ,K, For notational convenience we define 

C 4 diag(C^,...,Cj^) 

R = dlag (JR 

Comments; 

The generalized conic sector error bound (18), though somewhat 

1 

abstract, is fairly easily related to meaningful quantities. 

Lemmas A4 and A5 in the Appendix relate condition (18) to simple : 

Euclidean norm bounds for multi-input- multi- output (MliViO) NTV and LTI 
For example, if is a stable SISO LTI element, then it follows from 
Lemma A4 that the simple frequency domain condition 


(19) 

( 20 ) 
( 21 ) 
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1 - C|, (jw) I ^ < 1 r^(ju)) I ^ 


(22) 


- e 


for some e > 0 and all u), (see Fig. 2b) implies that (18) is satisfied for any 
stable SISO L"’I nonanticipative S and R satisfying 


c.(jm) 

= c.(joj) 

(23) 

R.(» 

Sj,(j“>) 

= r.(ju)) . 

(24) 


If H . is a SISO memoryless NTV element defined by 

M X 

y^(t) = h.(x.(t),t) (25) 


and if for sor. ; 5 £ and some e > 0 


t) - 
a 


2 


< 



V a 0 , 


(26) 


then from Lemma A5 it follows that (18) is satisfied for any constant 
C^(s), R£(s), S^(s) satisfying (23) - (24); this is the usual sector nonlinearity 
condition (e. g, , [5-8])--see Fig, 2a. In more complicated situations where 
a component H* consists of an interconnection of several LTI elements 

X 

and/or memoryless NTV elements, the result of [21] frequently may be 
invoked to determine suitable G., R., andS.. 

~i ~ i ~ i 
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(a) Nonliviear component satisfying (Z6) 


(b) Nyqliist locus of LTI component satisfying (22) 
Fig. 2 SISO Components 
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IV. MAIN RESULT 


it 


Our main result is now stated. 


Theorem I (Multiloop Circle Criterion) 

I 

Let C, R, and S have respective proper rational transfer 
function matrices G(s), C(s), R(s), and S(s); let R~^ and S' ^ exist; and 

ft* M 

let C, R, R"^, S, andS'^ be nonanticipative and L»^-stable; let H . 

(i = 1,...,N) satisfy the condition (18). Suppose the feedback system (11) - 
(12) is nonanticipative and L^^- stable in the special case where 
H . = C. V*i = 1, . . . /N. Then, a sufficient condition for the system 
(11) - (12) to be L_ -stable for every collection H. (i = 1, , . . ,N) satisfying 

w 6 1 

(18) is 


G(jw) (l + C(jio) G(juj))"^S^(ju>)) < 1 Viu. (27) 

Further, when G(s) is a square matrix and is invertible almost everywhere 
on the jo)-axis, then the above condition can be expressed as 


min 


(S(j(o) (C(ju)) G"^(jo))) R'^juj)) > 1 


(28) 


for all u) at which G"^(ju)) exists. 

Proof: From Lemma A3 and (27), it follows that 


□ 


H ^ diag(Hj, 


Hj^) strictly inside L^^ - Cone(G, R, S) , (29) 


From Lemma A2, it follows that (27) and (28) are equivalent when G" ^ (jiu) 
exists almost everywhere and that 
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(30) 


(-G)^ outside L, -Cone(C, R, S) , 

L^g-stability of (1X)-(12) follows from the conic sector stability theorem, 

Theorem A4, ^ 

□ 

Remarks 

It is also possible to prove Theorem 1 by applying ParsevalVs 
theorem to verify that the conditions of the well-known "small gain theorem" 
(e. g. , [8]) are satisfied by the transformed system defined by Eqns, (A5)- 
(AlO) oi the Appendix, We consider the present proof more appealing 
because it stresses the direct connection between the conditions of 
Theorem 1 and the simple conic sector conditions of Theorem A4, 
just as Zames' proof in [22] of the well-known circle criterion stresses 
the direct connection between circle theorem conditions and the simple, 
but less general, conic sector conditions of [5, Theorem 2a], 


V. DISCUSSION 

There are essentially two main conditions which miist be satisfied 
to conclude stability from Theorem 1: (i) The system must be stable when 
the uncertain components are replaced by the respective LTI approxi- 
mations C.; and, (ii) the frequency-domain condition (27}{or (28))must be 
satisfied. The former condition can be verified a variety of ways; for 
example, one may check that the roots of the characteristic equation 

det (I -f C(s) G(s)) = 0 (31) 
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4 

all have negative real parts; alternatively, one may apply the multivariable 

Nyquist criterion, checking that the polar plot of the locus of dot (C(ju)) G(juj)) 

encircles the point -1 + jO exactly once counterclockwise for each unstable 

open-loop pole of C(s) G(s) (multiplicities counted) f3]. The latter condition 

(27) (or (28)) requires that one plot the variable {• ) (or (♦ )) verses 

max min 

uj and verify that the appropriate inequality holds for all u), 

In the special case in which there is a single SISO nonlinearity 
h(x, t) both of the conditions of Theorem 1 can be verified by inspection of 
the polar plot of G(jU)) V'UJ. In this case the conditions of Theorem 1 
become precisely the conditions of the well-known circle stability criterion 
(cf. [7]“[8]). It is this which motivates us to refer to Theorem 1 as a 
"multiloop circle stability criterion"--despite the fact that in general no 
circles are employed in verifying its conditions. 

One can b.terpret tb,e uncertainty bounds (R^, S^) as specifications 
for the gain margins and phase margins of tiie system (11)- (12). For 
example, if the li . are LTI and SISO and if H.(s) = C.(s) (i = 1, . . . , N), 
then it follows from (22) that under the conditions cf Theorem 1 the system 
will remain stable despite variations in the individual component gains of 
magnitudes great as 1 r^(ju)) j = j Rj^(juj)/S^(ju)) ] , even when the variations 
occur Simultaneously in all components. So, for example, the system can 
tolerate simultaneous gain variations or phase variations of at least 



inf 

OJ 


20 log 


r.(jiu) 


db 


( 32 ) 


^If C(s) G(s) has any "decoupling zeroes" (i. e. , uncontrollable or 
unobservable poles, then these will not be roots of (31) and one must 
check separately that these poles have negative real parts- -cf, [3], 
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or 




a res in 


2C.(juj) 


{33) 


in each o£ the respective component feedback loops; i, e. , the system has 
gain margins of at least and phase margins of at least 0j^, at the inputs 
to the respective components Cj^(s) {i = 1, . . . ,N). The quantity 


is the amount by which the uncertainty bounding matrices R^(jiu) can be 
simultaneously increased without violating the stability conditions of 
Theorem l“”k^{ju)} can be viewed as a lower bound on the amount by which 
the system (11)-(12) exceeds the stability margin specifications (18) at 
tiach frequency ^ . 

In general, the stability conditions of Theorem l--and the estimate 

(34) of excess stability margin k„ — will be conservative. This conserva- 

m 

tiveness can usually be reduced by substituting weighted uncertainty 
bounding matrices R^(s), S^(s)) for the original matrices 

(R(s), Sj^(s)), Further, if it happens that is linear time-invariant for 
some !, then as a consequence of Lemma A4, the corresponding cv. 
may be replaced by a frequency-dependent a*(jai) (provided that for some 
k<<» and some e>0, k>|o^(jw)| >s for allgu). Iterative numerical methods 
would be required to enable one to efficie ntly compute the "optimal" 
weightings (i. e, , the weightings leading to the least conservative stability 
conclusions) . We hasten to add that the Idea of using constant weightings 
to reduce ronservativeness in multilopp nonlinear input- output stability 
results is not new; M-v.iatrix tests provide a simple but conservative 
method to implicitly ensure the existence of constant 
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weightings (see [I’O], [ll]); other results have been stated in which constant 
weightings appear explicitly (e, g, , Moylan and Hill [12, Thms. 5 and 6]) . 

The results of [12] involving explicit weightings may be viewed as a 
special case in which the matrices (C{ju)}, R(ju>), S(ju))), the interconnection 
matrix G(jw), and the weightings of. (jw) are not permitted to be frequency- 
dependent, We emphasize that the advantages offered by frequency- 
dependent C(j(u), R(jis), S(juj) and a^(juj) are crucial in stability margin 
analysis where it usually is necessary to be able to characterize tolerance 
of dynamical modeling errors of frequency- dependent magnitude. Allowing 
the matrix G(juj) to be frequency dependent eliminates the need for 
incorporating the dynamics of G in additional dynamical H.'s, thereby 
reducing the dimension N of (G, C, R, S); this in general leads to less 

M ^ 

conservative stability conclusions from Theorem 1 and also broadens its 

scope of applicability (since G need not be stable under the conditions 

of Theorem 1 whereas each must). 

We note that Theorem 1 is fairly broad in its scope of applicability. 

The transfer matrix G(s) may be non-square and need not be open- loop 

stable (though the H. must be). The H. operators may be raulti-input- 

~ 1 — 1 

multi-output and need not have equal numbers of inputs and outputs. Unlike 

some previous interconnected system results, no condition is imposed 

2 2 2 2 

requiring either c^^ - r^ > 0 Yi or c^ - r^ <0 Yi. The operators 
(C., R., S .) defining the conic sector condition (18) may be dynamical. 
Further, with the aid of Lemmas A4 and A5 and the aid of results such as 
in [21], it is practical to verify the conic sector condition (18), even for 
multi- input- multi- output dynamical nonlinearities. 
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VI, CONCLUSIONS 


With a view towards developing a stability criterion well suited 
to the problem of multiloop feedback stability margin analysis, nonlinear 
input-output stability techniques generalizing the circle criterion have been 
re-examined. The s tabiHty margin implications of existing results have 
been stressed and an improved result has been generated allowing one to 
take account of the frequency-dependence of the magnitude of system 
modeling errors which commonly occurs in situations involving imprecisely 
modeled dynamics- -e. g. , singular perturbations, hysteresis, etc. 

Theorem 1 together with the related Lemmas A1 and A2 provide verifiable 
sufficient conditions for the stability of multiloop feedback systems using 
only crude conic sector bounds on system parameters, subsystem 
frequency responses, and nonlinearities. Potential applications include 
the testing of system integrity in the presence of actuator and/or sensor 
failures {cf. [2]) and the characterization of frequency-dependent gain and 
phase margins for multiloop feedback designs subject to multiple singular 
perturbations and dynamical nonlinearities leading to simultansous 
frequency-dependent variations in gains and phases in the feedback loops. 
The main result. Theorem 1, also plays a key role in a result 
described in [21] for generating conic sector bounds to characterize the 
sensitivity of multilocp systems to large dynamical modeling errors of 
frequency-dependent magnitudes in a manner similar to the wa/ classical 
M-circle and Nichols diagram tecliniques enable one to quantitatively 
gauge the effect of open-loop gain variations in single-loop, unity-feedback 
systems. The result of [21] also can be useful in determining the conic 
sector bounds (C, R, S) required by Theorem 1. 



VII, APPENDIX 


In this appendix several results are stated which are needed in 
connection with Theorem 1. Theorem A1 and Lemmas A2 and A3 are used 
in the proof of Theorem 1. Lemmas A4 and AS are useful in verifying the 
conic sector conditions (18) for memoryless nonlinear H, and for linear 
time-invariant dynamical H 

1 

We note that while results similar to Lemmas A2, A4, and A j have 
been presented in various forms elsewhere (e. g. , [6]-[8], [13], [19], [20]), 
the very general case considered here (admitting, for example, dynamical 
and multi-input- multi- output C, R, S) is new, as is the explicit appearance 
of the term c llytl^ in. the "strict" conicity condition (5)-(6). The differ- 
ences are sufficient to mandate the inclusion here of proofs for these 
Lemmas, 

Theorem Al (Conic Sector Stability Theorem) 

Consider the feedback system 


H X 

(Al) 

-G(y + v) + u 

(A2) 


where x, y, u,veL 2 ^ and G, H: L^^ — L^^j (u, v) is the "input" and (x, y) 
is the "output". If L^^-stable linear operators C, R, and S can Le found 
such that 

H strictly inside L 2 g-Cone (C, R, S) (A3) 

and 

(-G)^ outside L^^-Cone (C, R, S) , (A4) 
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then the feedback system (Al)-(A2) is L< 2 q- stable. 

Proof; This result is a special case of the "sector stability theorem" 

([13, p. 65], [19, Thm. 6.1]). 

P.emark 

Theorem A4 also can be proved by applying the small gain theorem 


(e. g. , [8]) to the "transformed" system 

y = S (H - C) R“ ^ X + V (A5) 

X = -R G (I + C G)"^ y + ti (A6) 

whei'e 

y = S (y - C(x-u) + v) (A7) 

X = R X (A8) 

u = R u (A9) 

V = S (v + C u) , (AlO) 


provided S"*^, G(I + C G)"\ and R“ ^ exist and are L- -stable and non- 

anticlpative. It can be shown that under the conditions of Theorem A4, the 
two operators S (H - C) R ^ and RG(I+CG) ^ each have 
less than one. Thus, Theorem A4 may be viewed as a characterize.tion 
of the improvements on the small gain theorem obtainable by use of the 
nonlinear dynamical transformations (A5)-(A10); such transformations 
have been described as "loop- shifting" and "multiplier" transformations 
(e. g. , [8]). 

Lemma A2 (LTI Outside Conicity) 

Let G, C, R, S be linear time -invariant operators with respective 



proper rational transfer functions G(s), C(s), R(s), S(s)» Suppose that 
S“^(s) exists anqi has a proper rational transfer function matrix with no 
poles in Re{s) > 0, Suppose that R, G(I + C G)”^ and S"^ are stable 
and non-anticipative. Then, 

(-C)^ outside L-^-Cone(C, R, S) (All) 

if and only if the following condition holds for all real ou 

G{juj) (l + C(ju)) G(jw))'^ S’^jw)) < 1 (A12) 

almost everyv/here. 

When G“^(s) and R"^(s) exist almost everywhere on the jiu-axis, then 
condition (A 12) is equivalent to 

almost everywhere, 

Proof; 

It is trivial to see that (A12) and (A13) are equivalent when G“^(ja>) 

“ 1 

and R (j«J) exist, since for any invertible matrix A 

a . (A“^) = 1/a (A) . 

mm ' ’ max' ' 

Suppose that (A12) holds. Let (x, y) by any input-output pair 
satisfying X = - Gy; let 

y = S (y - C x) (A 14) 
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and let 



y(t) , if 0< t < T 
0 , otherwise 


(A15) 


Let denote the Fourier transform of Note that from (A12) it 

follows that for all 

11 - llR(ju>) G(j<») (I + c(ju)) G(juj))-^ S'l(j<B) 

Now, 

Wlx-Wr = + 11 ^ 

by the nonanticipativeness of n 
R, G ( I + C G)"\ S"^ 

CO ^ 

< r 11 R G (I + c G)"^ S“^ Y_) (t) 11 ^ dt 
0 

:r ^ J 11 R{j(u) G(ju)) (1 + C(jU)) G{juj))"Vs"^joj) Y^(ju>)f dU) 

. CO 

^ by Parseval's Theorem and the hypotheses that ^ 

R, G(I + C G)"^ and S" ^ are L,^-stable. 
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< -^ J 11 Y^(juj) 11^ duj 




- 11 11t = I! y 11t 


= il S(y - Cx) 11 


Z 


(A17) 


which proves (All) is implied by (A 12). 

Conversely, suppose that (All) holds. Let Yq and be arbitrary. 
Consider the L^^-stable feedback system 


X = -Gy 


-1 


y = C X + S y 


(A 18) 
(A19) 


Let y(t) -*• Yq e . Then (letting T •♦ « in (7)), it follows from (All) that 


!' R(ju)) C(Ju.) (I + C(jc») G(i..')>‘‘ Y|j 111 


S 11^0 111 


(A20) 


and hence (A 12) holds. 

Lemma A3 (Composite Operator Conicity) 
Let 

H = diag (H , Hj^) 

C = diag (C , pjvj) 


□ 


(A21) 

(A22) 
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R s diag (R . . . ,Rj^) 


(AE3) 


S = diag (S . . . »Sj^) . 


(A24) 


If for all i s 1 , . , . , N 


H. strictly inside L^^-Cone (C., R S , (A25) 

then 

H strictly inside L- -Cone (C, R, S) , (A26) 

Proof; This is a special case of the results in [13, p, 70] and [19, 

Lemma 6.2 (vi)], ^ 


Lemma A4 (LTI Conicity) 

Let H, C, R, S be nonanticipative L 2 ^“Stable linear- time- invariant 
operators with respective rational transfer function matrices H{s), C(s), 
R(s), S(s). Suppose that R (s) exists and has no poles in Re(s) > 0. 

Then 

H strictly inside L^.-Cone (C, R, S) (A27) 

pm (U C mp pm pm 


if and only if 


!1 s(ju)) (H(juj) - C{jm})x(j(s)'.l^ < 


2 


!l X{jio) 11 

E 


- e 


i! x{ju)) 1 


2 

E 


(A28) 


for all X(juj), all a) , and some s > 0* 
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« 1 

Proof: Let R, denote the stable nonanticipative LTl operator having 

_ 1 

tranefer function matrix R (s). Suppose that (A28) holds and let 


- 


(Rx)(t) , if t< T 


(A29) 




, if V. > T 


and letx(jtu) denote the Fourier transform of x_{t). Then, for all y = Hx 
j T ~ 

we have 


cx) llj = |l_s(Hx- cx)|l^ 


II S(H - 9 == II T 

(by linearity) 

II S(H - C) R“^Rx H 


(since R” ^ exists) 


- 1 II ^ 


S(H - C) R“^ ! 


(by nonanticipativeiiess) 


< J II (S(H - C) R"^ x^)(t)l| dt 
o''" ^ 

-1 


(the integral exists since S, H, C, and R are stable) 




i-) • . 


w - 

= J 11 S(jw)(H(juj) - C(ju)))R-^(joi) 


X^(jw) IP d(JJ 


E 


(by Parseval's Theorem) 
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if- ** 


•» i-/ ♦'•i •‘f •* 


< ilx^(juj)l! - E |lR~^(ja)) x^(jw) P dw 

S E 


(By Ar^)) 


CO 2 

“ IllLWll dt- « f ll(R-'xJt) ll^dt 

0 ^ E J ~ ^ E 


(by Parseval's Theorem) 


lix (t)li^dt- e j’'ll(R-^~x )(t)!l^ dt 
Jq T E Jq - " E 


11 ^ - E il X^ 11 ^ 


R xli - £ 11 llx 


(by nonanticipativeneas of R 


-1 


= II Rx II, 


e x 


< R X 

M * 


e'(|UL^+ IlHxL^ ) . 


where 


and 


e = 


1 + a 


Hx 


a - 


sup 
X, T 


■T 


llx 11 , 


< o 


(since H is stable). 


(A30) 


(A31) 


(A3 2) 
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Conversely suppose (A27) holds. Let and u>q be arbitrary. 

Then, letting. x(t), -* Xq cos {(u^t + 56) and T it follows that 


ti 

11 SOalg) - COul;,)) X|J II 

2 2 2 

< 11 R(ju.g) 11 - e ( II II + II H(ju.g) X„ || ) 

£) £j L 

< |lR{jmg) X^l|2 - ^ l|X<j ll^ . 

£ £ 


Lemma AS (Nonlinear Time- Varying Conicity) 

Let h(x(t), t) be any function of x(t) and t and let H be given by 

(Hx)(t) = h(x(t), t) . 

A/ 

Let C, R, and S be constant matrices and let C, R, S be the operators 
defined by 


(Cx)(t) = C x(t) 

¥-x 

(Rx)(t) s R x(t) 

-¥-x 

(S y)(t) = S y(t) . 

^y . 


Suppose S ^ exists, then 

. I • 

H strictly inside L* -Cone (C, R, S) 
■»V \iy- .il ® ~ ~ ~ 

li 

if and only if ' ' 

- Z8- 


(A33) 

□ 


(A34) 


(A35) 
(A3 6) 
(A37) 


(A3 8) 



V x(t) . (A39) 


Z Z 

11 S(h(x(t),t) - Cx(t)) 11 < llRx(t)ll 

E E 


2 


I! ^(t) 11 

E 


□ 


Proof: Let y(t) = h(x(t), t). 

Supf^ se (A3 9) holds. Then, 


11 y(t) 11 < a \\ x(t) II 


(A40) 


where 


‘^max 

" = a ' ": W 

min ' 


(A41) 


Thus , taking 


e' = 


1 + a 


(A42) 


we have that 


S (y - C x) 


= [ 11 S(h(x(t),t) - Cx(t) 11^ dt 

„ E 


r.T 


R x(t) 11 - f. 11 x(t) II dt 

0 E E 


J 11 11 - e'(llx(t) 11 g + 

0 ^ 


y(t) 11 dt 
E 


- t ll yll^^) . 


(A43) 
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Conversely, when (A38) holds, then taking x(t) to be the constant 
function x(t) = .tCq we have that for some e > 0 

llS(h(x„.t) -Cx^) f 

= ^ ■ II 5 f?’'- IIt 

< Kii?-iix-ii-'=iO 

= II RXg 11^ - e 11 X(, 11^ . 


(A44) 

□ 
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